Confinement, Chiral Symmetry Breaking and complex mass in QED2+1 in Minkowski 

and Euclidean spaces 



Vladimir Sauli 1,2 and Zoltan Batiz 1 

'CFTP and Dept. of Phys., 1ST, Av. Rovisco Pats, 1049-001 Lisbon, Portugal 
2 Dept. of Theor. Phys., INP, Rez near Prague, AVCR 

Without any analytical assumption we solve the ladder QED2+1 in Minkowski space. Further, we 
transform Greens functions to the Temporal Euclidean space, wherein we show that in the special 
case of ladder QED2+1 the solution is fully equivalent to the Minkowski one. The obtained complex 
fermion propagator exhibits confinement because of the lack of the pole at the real timelike q 2 axis. 
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I. INTRODUCTION 
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Quark confinement in Quantum Chromodynamics (QCD) is a phenomenon of current interest. Due to the fact 
that QCD is not easily tractable, the toy models which exhibit QCD low energy phenomena -the confinement and 
chiral symmetry breaking- are often investigated. In this respect 2+ld Quantum Electrodynamics (QED2+1) has 
these similarities with QCD in the usual 3+ Id Minkowski space. Based on the Euclidean space study of QED3, the 
chiral symmetry breaking for a small number flavors has been proposed for the first time in [l[. Since the scale of 
dynamical chiral symmetry breaking, being characterized by a fermion mass in the infrared -M(0)-, is one order of 
• magnitude smaller then the topological dimensioned coupling e 2 , the Schwinger-Dyson equations (SDEs) provide a 
Oh 1 unique powerful framework for the nonperturbative study, see e.g. most recent studies [IJl]. The importance of the 
unquenching effect in QED3 for an increasing Nj has been recognized a long time ago and reinvestigated in 

SDEs framework lately 0, H 0> HI • Particularly, confinement in relation with dynamical complex pole generation in 
. fermion propagator has been discussed in 0] ■ Further analogy with QCD in finite temperature and chemical potential 



has been explored 



It is noteworthy that QED3 is of current interest as an effective theory of high-temperature 



£T) . superconductors U TT|, E3, E3] , Mott insulator and graphene [HI, EH 

■ However, the all aforementioned studies have been done in the standard Euclidean space. That is, after performing 
the standard Wick rotation [T7| of the timelike components of the momentum variables (internal integral momentum 
as well as external one, explicitly pf = —ip^ 1 , the measure io9p M — —d 3 p E ). It is assumed and widely believed, that 
the Green's functions for timelike arguments can be obtained after analytical continuation of the functions calculated 
in Euclidean space, where it is supposed the Euclidean solution itself represent the correct Minkowski solution for the 
spacelike arguments (for some discussion see also Sections 2.3 and 6.3 in 18]). Therefore, to shed a new light and for 
the first time, we solve fermion SDE directly in 2+1 Minkowski space. The so called ladder approximation of electron 
SDE is introduced in the the Section II. 

In the Section III, by using hyperbolic coordinates, we solve the SDE directly in the Minkowski space. While for 
timclike p 2 it provides selfconsistent gap equation, i.e. the integral equation with the mass function of the same 
arguments in both sides of Eq., whilst for spacelike Minkowski subspace the solution is naturally constructed from 
necessarily known timelike solution. Furthermore, we derive the SDE in the Temporal Euclidean (Et) space in the 
Section IV and prove that this exactly leads to the original Minkowski formulation for timelike momenta. Recall here, 
Et space metric is obtained from Minkowski one by the multidimensional Wick rotations, but now instead for the 
time component, it is made for all the space coordinates of the Lorentz three vector [jjjj . 

The numerical results for various ratio of the coupling and the electron mass are presented in Section V. The 
imaginary part of the mass function is automatically generated for a coupling strong enough and resulting fermion 
propagator has no mass singularity in the timelike region. Such electron electron can never be on-shell and thus 
never observed as a free particle, in other words it is permanently confined @, H9, [H, [Hj]. We further discuss the 
(in-)validity of standard Wick rotation and conclude in the Section VI. 



II. FERMION SDE IN QED2+1, QED3 IN LADDER APPROXIMATION 



In our study we employ Minkowski metric g^ v = diag{\, —1,-1), in order to properly describe chiral symmetry, we 
use the standard four dimensional Dirac matrices such that they anticomutation relation is {7^,7^} = 2g M „. With 
these conventions the inverse of the full fermion propagator reads 



S 1 (p) = ^-m-S(p) 
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—-,G^(k-p)r fl (k,p)S(k) 7v (2.1) 

We consider the explicit chiral symmetry breaking mass term ol the form mifiip so parity is conserved. In this case 
the dressed fermion propagator can be parametrized by two scalar function like 

S(P) = S V (P)^S S{P )= M{p) 1 _ B(p) . (2.2) 

The full photon propagator G and the electron-positron-photon vertex T satisfy their own SDEs (for their general 
forms see |18j ) 

The ladder approximation is the simplest selfconsistent approximation which approximate the unknown Greens 
functions be their free counterpartners, i.e. T mu — j mu and the photon propagator in linear covariant gauges is 

n _ + (1 ~~ 0~W^ , v 

m v p • \ z - 6 ) 



III. DIRECT MINKOWSKI SPACE CALCULATION 



In general QFT the Greens functions are not real functions but complex tempered distributions. In perturbation 
theory these are just real poles (together with its Feynman ie prescription) of the propagators, which when coincide 
in the loop integrals, produce branch cut starting at the usual production threshold. At one scalar loop level, the 
two propagators make the selfenergy complex above the point p 2 — (Mi + M2) 2 , wherein Mi, M2 are the real masses- 
in fact the positions of these poles. Depending on the masses of the interacting fields, the real pole persists when 
situated below the threshold or we get non-zero width and the free particle becomes resonance with finite lifetime. 

In strong coupling quantum field theory the mechanism of complexification can be very different (however the 
mixing of both mechanisms is not excluded). Here we simply assume that there is no zero at the inverse of propagator 
for real p 2 , thus ie factor is not necessary and we integrate over the hyperbolic angles of Minkowski space directly. For 
this purpose we have to consider the propagator function as the complex one for all p 2 . A convenient parametrization 
of the complex fermion propagator functions S s and S v can be written as 



S s (x) 



B(k) 



S v (k) 



A 2 (k)k 2 - B 2 (k) 
= Rb [Cg| - T 2 A )k 2 ~R 2 B - £|] + 2R a T b Ta k 2 

D 

■ ?b [{R\ - T 2 A )k 2 +R 2 B + r|] - 2R B R A T A k 2 

+ * — J) ' ( 3 ) 

A(k) 



A 2 (k)k 2 - B 2 (k) 
R A [(R 2 A + T 2 A )k 2 - + £|] - 2R B T A T B 
D 

. T A [-(R 2 A + Y 2 A )k 2 -R 2 B + r|] + 2RaRbT b r _ 

+ 1 Y) ' \ 6 - z ) 

where R a ,Rb (Fa^b) are the real (imaginary) parts of the functions A, B and the denominator D reads 

D = {[R 2 A - T\]k 2 - [R% - r|]) 2 + A(T A R A - T B B) 2 . (3.3) 

In order to be able to compare between Minkowski and lately considered Euclidean spaces, the transformation in 
use should leave the Minkowski spacetime interval 

s = t 2 - x 2 - y 2 (3.4) 

manifestly apparent (i.e. s could be a variable of the integral SDEs). 

To achieve this we will use 2+1 dimensional pseudospherical (hyperbolic) transformation of Cartesian Minkowski 
coordinates. The obstacles followed by Minkowski hyperbolic angle integrals when going beyond A — 1 approximation 
restrict us to the Landau gauge wherein the A = 1 is the exact result in Temporal Euclidean space. 
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In momentum space our convenient choice of the substitution is the following: 



/ 



d 3 kK{k,p) = 



drr 



cosh< 



d6 



da < sinha 



k a = —r cosh a 
k x = —r sinha sin# 
k,, = —r sinha cos 9 



sinh a 



k a = r cosh a 
k x — r sinha sin# 
k,., = r sinh a cos 9 



k Q = —r sinha 
k x = —r cosh a sin£ 
k y = —r cosh a cos( 



cosh ( 



k Q = r sinh a 
k x = r cosh a sin 9 
k y = r cosh a cos 9 



K(k,p). 



(3.5) 



Notice, the integral boundaries are universal for all the subregions of Minkowski space, the first line corresponds to 
the integration over the timelikc 2+lmomentum where we have 



k 2 =r 2 >0, 



k 2 — k 2 k 2 

where the left term corresponds to the negative energy interval ko < 



(3.6) 



k 2 and the right term corresponds to 



the positive k^ > + \/k 2 + ky. The second line stands for the spacelike regime of the integration 

k 2 = -r 2 < , 



(3.7) 



wherein the left term corresponds to the energy component interval fco = {—\ k 2 + k 2 0), while the right term in the 



second line stands for positive k Q = (0, yk 2 + k 2 ) subspace of the full 2+1 dimensional Minkowski space. Functions 

V in Rel. (|3.5p represents the integrand of SDE. 

The functions A, B are Lorentz scalars, thus they can depend on p 2 only. We freely take the simple choice of 
timelike external momenta as p^ — (p, 0, 0) which leads to the following a integrals for the timelike part of internal 
momenta: 



da- 



sinh ( 



P 



2pr cosh a 



sinha 1 
da— — = — In 

r z _|_ pz _ 2p r cosn a pr 



r — p 



r +p 



(3.8) 



The contribution stemming from spacelike part of loop momenta gives zero because of negative and positive time 
volume contributions, although each infinite separately, they cancel each other. Note, the 9 integrals contributes by 
simple 2tt prefactor. 

Integrating over the angles we can see that at the level of our approximation, the SDE separate for spacelike and 
timelike regime of the threemomenta. For timelike p we get for the function B 
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(3.9) 
is derived 



where £ is a gauge parameter (and where we trivialy changed r of (|3.8p to k) Stressed here, the Eq 
without any requirement of analyticity for the the propagator and the kernel. 

For the external spacelike Lorentz threevector of momenta the a integration over the spacelike regime gives zero. 
Taking for instance p^ — (0,p, 0) this can be most easily seen by the inspection of the integrals 



da- 

a 
b 



coshf 



V + p 
da 



2pr cosh a sin t 
1 



da- 



sinh ( 



2pr cosh a sin 9 



1 



b cosh < 



b cosh < 



(3.10) 



where a 



, b = 2pr sin 9. The integrals in bracket (|3.10p can be evaluated by using of the following formula: 



da- 



1 



1 



da- 



b cosh a 
1 



y/a 2 + b 2 
1 



In 



; ln 



a - 


h& 


f Va 2 


-b 2 


a - 




-Va 2 


-b 2 


a 4 


-b- 


f Va 2 


+ b 2 



b sinh a Va 2 + b 2 a + b- yfa 2 + b 2 



(3.11) 
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The remaining contributions stem from the combination of the external spacelike and the internal spacelike mo- 
menta, putting now a = p 2 + r 2 , b = 2pr sin 9 we can get 

1 sinh a f sinh a 

r 2 + p 2 + 2pr sinh a sin 9 J r 2 + p 2 — 2pr sinh a sin 9 

roc / l j x 2 In (V 1 + (&/a) 2 - tya) 

= -r / rfc M r-r-; , . = = -: J - 3.12 

6 J V a + ° smn a a ~ ° smn a / b sin 9 a/1 + (6/a) 2 

where the definite integrals in the second line are given by second line in Eq. (|3.11|) . 

Adding the all together we get for p 2 < (for £ = and A = 1) the following integral expression for the function 

B: 



2jv 



ln(VTT^2-z) B(r 2 ) 



B ip )=rn + i 7 -^ I dr\ d9-^ ,„ ^ 7 , „,,,.,, 



(27t) 3 7 J sin0\/T+~? r 2 -B 2 (r- 

2 \] —p 2 r sin 

2 i 2 ' 

r + p 



(3.13) 



To conclude, we stress the main difference when compared to the standard treatment. Here this is the timelike part 
of Minkowski subspace where the results are most naturally obtained. Quite opposite to the standard approach where 
Minkowski solution is constructed by continuation of the Euclidean result, here the solution for the spacelike argument 
is non-trivialy made from the timelike soluiton which must found as a first. Actually, the solution at timelike domain 
given by Eq. (|33|) is needed in the rhs. of (|3~T3)) . 

The validity of standard Wick rotation is highly speculative topics in the literature. Here we know the solution 
directly in Minkowski space and the comparison with the Euclidean solution is straightforward and easy task. Assum- 
ing the Wick rotation is valid, we could get the Minkowski solution equal to the Euclidean one at spacelike domain of 
momenta. We anticipate here and it will be explicitly shown in the next Section that the mass function B is actually 
complex for all real timelike p 2 at large window of studied parameters m and e. Hence the Minkowski solution B(0) 
being complex, is not the one literally known from the Euclidean studies, where -Be(O) is an always purely real num- 
ber. These solutions do not coincide in the Minkowski light-cone (Euclidean zero) and in other words: commonly used 
strategy based on the analytical continuation of the Euclidean (spacelike) solution to the timelike axis is wrong and 
should be abandoned in the case of QED2+1 theory. In our simple model, the "continuation" is not the analytical 
continuation of necessarily holomorphic function, but is given by Eq. (|3 . 1 3[) . This should be reconsidered in any other 
specific model -approximation and truncation of SDEs- system. 

We do not know whether there is a unique general rule to continue between spaces (we mean a simple prescription, 
which do not require the solution of SDE for completely general complex arguments). However, before presenting 
the numerical solution of the Eq. (|3.9[) . we will exhibit the equivalence of QED3 defined in the ET space, where the 
metric simplifies and allows to consider more complicated approximations of SDE. 



IV. PROOF OF AN EQUIVALENCE OF QED FERMION SDE FORMULATED IN MINKOWSKI AND 

TEMPORAL EUCLIDEAN SPACE. 



In this Section we will show that the Temporal Euclidean formulation of our considered SDE exactly agrees with 
the one derived in Minkowski space. 

In the paper 19] it was proposed that N-dimensional analog of Wick rotations performed for space components 
of Minkowski N + 1-vector can be partially useful for the study of a strong coupling quantum field theory. In 
even dimensional 3+1QCD it is the nonperturbative mechanism responsible for complex mass generation which is 
responsible for the absence of real pole type singularities in Greens functions evaluated at real their arguments which 
thus makes the nonperturbative calculations feasible there. 

In odd dimensional theory, like QED2+1 we study here, the complexification of masses and couplings can be quite 
naturally expected because presence of i in the measures of the integrals defined in ET space. To see this explicitly, 
let us consider the momenta as complex variables and let us assume that there are no singularities in the second 
and the fourth quadrants of complex planes of k x , k y in the momentum integrals (see I2.ip . Deforming the contour 
appropriately then the aforementioned-mentioned generalized Wick rotation gives the following prescription for the 
momentum measure: 

k x , y — > ifc 2 , 3 , 

i I d 3 k^-i I d 3 k Er , (4.1) 
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which, contrary to our standard 3 + 1 space-time, leaves the additional i in front. 

In Et space the singularity of the free propagator remains, for instance the free propagator of scalar particle is 



1 



• is 



with a positive square of the three-momenta 



■Pi 



pi-. 



(4.2) 



(4.3) 



thus formulation of the weak coupling (perturbation) theory, albeit possible, would not be more helpful then the 
standard approach (Wick rotation). 

The advantage of the transformation to Et space becomes manifest, since the fixed square Minkowski momentum 
p 2 = const hyperboloid with infinite surface is transformed into the finite 3dim-sphere in Et space. The Cartesian 
variables are related to the spherical coordinates as usually: 



&3 = k cos 8 
k± = k sin 9 cos < 
fc 2 = k sin 9 sin (j 



(4.4) 



Making the aforementioned 2d Wick rotation, taking the Dirac trace on E and integrating over the angles we get 
for the function B 
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~P 



S s (k) 



(4.5) 



which is the same equation we derived directly in Minkowski space (|3.9|) . 

Making the trace Tr /pY,/{Ap 2 ) and integrating over the angle we get for the renormalization wave function in Et 
space 



m 
i 



i 



i 



4tt 2 

2 i 



k 2 

dk Sv (k 2 )[~i + {\-i)i] 



k 2 



2pl 



In 



P 



k+ p 



(4.6) 



with the propagator function S v defined by (|3.2[) . The first term in the bracket [] follows from the metric tensor while 
the second one proportional to (1 — £) stems from the longitudinal part of gauge propagator. We can see that like in 
the standard ES formulation we get A = 1 exactly in quenched rainbow approximation in Landau gauge £ = 0. 



V. NUMERICAL SOLUTIONS FOR TIMELIKE p 2 

Before describing our result, at this place it is noteworthy to mention the study of complex singularities which 
were found when the momentum was continued to the complex plane. The author of |7fl rotated integration contour and 
solve QED2+1 SDEs system for complex variables selfconsistently (i.e. the same ray pe 1 ^ with constant phase <fi was 
considered for the external and the internal variables). The location of the singularity was obtained by extrapolation 
of the complex solution, since at near vicinity of complex pole the numerical procedure was unstable. Recall here, the 
phase 7r/2 would transform standard Euclidean SDE to the Temporal space used here, however this solution was never 
achieved because of singularity found. As will be shown, the solution become stable again there, since we are enough 
far from the singularity. Actually, we have found the equation stable even when the pole is expected near real the 
timelike axis, since in our case the stability is guaranteed by a non-zero value of current fermion mass m in this case. 
Although we are not really interested in the location of the complex plane we can refer the result for the case e 2 >> m 
(only m = case were studied for all the all truncations of SDEs system in Q). The observed location of singularity is 
characterized by <j> ~ ir/4 and \p\ ~ M(0) where M(0) is the dynamical infrared mass. Apparently, the author did not 
continue the solution further to the timelike axis, because of crossing the singularity, however as we opposed above, 
doing opposite could be the correct procedure. One should start from the Temporal Euclidean solution, which being 
identified with the physical timelike Minkowski solution, the continuation to the spacelike regime could be performed 
afterwords (including the residua of the complex pole found). 

In this paper we present the solution of the mass function M — B in the ladder approximation in Landau gauge 
in (|4.5|) space, we do not perform the continuation to the spacelike axis, which remains to be done in the future. 
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We will consider the nonzero Lagrangian mass m and interaction strength characterized by charge e. QED2+1 is a 
superrenormalizable theory and as we have neglected the photon polarization it turns to be completely ultraviolet 
finite and no renormalization is required at all. We assume that the imaginary part of the mass function is dynamically 
generated and to get the numerical solution we split the SDE (|4.5[) to the coupled equations for the real and imaginary 
part of B and solve these two coupled integral equations simultaneously by the method of iterations. 

The complex mass generation observed herein is not only the simplest and most straightforward definition of quark 
confinement, but it has an advantageous technical consequence. The associated absence of the real singularity at p 2 axis 
ensures that, for all p 2 , the integrand is regular and hence the integral of SDE can be evaluated using straightforward 
Gaussian quadrature technique, i.e. there are no endpoint or pinch singularities on integration axis, except when the 
coupling constant is very small when compared to the bare mass m. In the later case the quantum corrections become 
negligible and the real branch point close to m is restored and heavy electron (or weakly interacting) QED2+1 is an 
unconfined theory again. 

To characterize complex mass dynamical generation it convenient to introduce dimcnsionless parameter 

e 2 

k=— . (5.1) 

2m v ' 

We set up the scale by taking m — 1 at any units. Pure dynamical chiral symmetry breaking (m = 0) is naturally 
achieved at large n limit. 

The phase <f>M of the complex mass function defined by M — \M\e l< ^ M is shown in Fig.l. for a various value of 
k. For very large n we get the dynamical chiral symmetry breaking in which case the obtained infrared phase is 
4>m{p 2 = 0, K = oo ) = 87.5° while more interestingly it vanishes for very small n. There is no confinement of fermion 
at all bellow some "critical" value of k , especially one can observe 4>m{p 2 — 0,k < 0.0191 ± 0.0001) = for one 
flavor ladder QED2+1. The absolute value of M is displayed in Fig. 2. for the various value of the coupling k. As k 
decreases the expected complex singularities gradually moves from complex plane to the real axis and the function 
develops expected threshold enhancement, however here it is much smaller when compared to QED3+1. 



VI. SUMMARY AND CONCLUSIONS 



We have presented the first analysis of the electron gap equation in Temporal Euclidean and Minkowski space. The 
main result, although based on the simple ladder approximation in given gauge, is the proof of the exact equivalence 
between the theories defined in Minkowski 2+1 and 3D Temporal Euclidean space. No similar is known about the 
standard Euclidean formulation and its relation with spacelike subspace of Minkowski space. 

The dynamical generation of imaginary part of the fermion mass can lead to the absence of real pole, providing 
simple scenario of confinement: there arc no freely moving fermions if k — e 2 /2m is not considerably smaller then 
unity. Minkowski QED2+1 has been shown to exhibit spontaneous chiral symmetry breaking -the mass function has 
nontrivial solution for m = and it is confining if the electrons are enough light. We expect robust quantiative 
changes when the polarization effect is taken account. 

The Temporal Euclidean space introduced recently and used here opens up a variety of questions. If there is 
any, what is the relation between standard Euclidean (e.g. lattice) formulation with the original Minkowski space 
confining theory? We clearly argue and explain in details in the Section III that the Wick rotation -the well known 
calculational trick in quantum theory- is based on an unjustified assumption for QED2+1. More interestingly, what 
are the outcomes when the same ideology is applied to the strong coupling even dimensional theories, e.g. to the most 
successful theory of the strong interaction in the nature: QCD [24|? 
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FIG. 2: Magnitude \M\ of the running mass M = \M\e l< ^ of electron living in 2+1 dimensions for different k, scale is m = 1. 



